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We study a two-dimensional effective orbital superexchange model derived for strongly correlated eg elec- 
trons coupled to f2g core spins in layered manganites. One finds that the ferromagnetic (FM) and anti- 
ferromagnetic (AF) correlations closely compete, and small changes of parameters can switch the type of 
, - magnetic order. For the same reason, spin order is easily destroyed with rising temperature, while altemat- 

' ^ ' ing orbital correlations can persist to temperatures where FM order has already melted. A scenario for the 

C ' AF phase observed in LaSrMn04 is presented. 

O I Published in phys. stat. sol. (b) 243, 277 (2006) 

Copyright line will be provided by the pubhsher 

> 

1 Introduction Manganites are characterized by a complex interplay of charge, spin, orbital and lattice 
Cp i degrees of freedom. Undoped and weakly doped single layer manganites Lai_j;Sri+j;Mn04 show antifer- 

^i I romagnetic (AF) order 1 1 1, in contrast to three-dimensional (3D) cubic Lai_2;Sra;Mn03 compounds |2|, 

where ferromagnetic (FM) ab planes are stacked antiferromagnetically in the undoped case (one electron 
C^ \ per site, x ~ 0). As for their 3D counterparts, the properties of the layered systems are strongly influenced 

by the orbital degrees of freedom of Cg electrons which couple to the spins and thus influence the magnetic 
C^ ' order. 

For a; = 0, the hopping of Cg electrons is blocked by large Coulomb interaction U, and charge fluc- 
, ^ tuations are quenched. They may be treated by second order perturbation theory which leads to virtual 

^ i djd'^ ^ d^d^ excitations by either Cg or t2g electrons, generating the spin-orbital superexchange 0. In 

Q I the present paper we intend to focus on the role played by the orbital degrees of freedom of Cg electrons 

O ' in stabilizing various types of magnetic order in the undoped single layer manganites, and the changes of 

spin and orbital correlations with increasing temperature. 

•T— I 
^> 

'V^ . 2 Effective spin-orbital model at finite temperature We study an effective spin-orbital superex- 

C^ I change model derived for the one-dimensional (ID) chain 0, generalized here to two-dimensional (2D) 

planes of layered manganites. The model depends on spin configuration S and takes the form, 

n{S) = Hj.+Hj + H,. (1) 

It consists of: {i) the superexchange Hji for the core spins formed by t2g electrons, (m) the orbital su- 
perexchange Hj for the Cg electrons, and {Hi) a crystal field term Hz- In the limit of large Hund's ex- 
change interaction, we restrict the Cg electron configurations to the subspace with their spins parallel to the 
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core spins at each site, and thus arrive at spinless fermions with an orbital flavor a = x, z, standing for 
\x) = -j=\x'^ — y^) and \z) = 4=|3z^ — r^) orbitals, respectively. Electrons with antiparallel spins are 
treated in second order perturbation theory in the same way as double occupancies. 

We start the discussion with the t2g core spin superexchange term Hji . In reality it depends on the 
total number of d electrons occupying two interacting Mn ions |5 1. However, in the undoped systems one 
finds d^ configuration at each site and one may replace this part of the superexchange by a Heisenberg 
Hamiltonian with a properly chosen exchange constant J' > favoring AF order: 

Hj,^J'Y,{S,-S,+i-S^). (2) 

i 

For convenience we use classical core spins Si of unit length, and compensate their physical value S = 3/2 
by a proper increase of J' . For classical spins Si parameterized by polar angles {i^i, (pi}, the spin product 
is given by (Si ■ Sj) — S^ (2\uij\'^ — l), where the spin orientation enters via 

Uij = cos(t9,/2) cos(i?j/2) + sin(79,/2) sin(t?j/2)e*('^^-'^*) = cos(%/2)e^X'^- , (3) 

depending on the angle dij between the spins and on the phase Xij^ which cancels at a; = 0. 

The superexchange term Hj describes eg orbital interactions, as derived recently for a ID chain J4l 
from the full coupled spin-orbital dynamics by replacing the spin scalar products by the actual averages 
Uij. The orbital superexchange Hj depends then on spin correlations via the factors {uy } in Eq. 0, 

^■' = l^ E (21^^^ I' + 3) {'^T^TJ - in,n,) - ^ J J^ (l " l"«.f )^<^K- (4) 

{ij}\\ab {ij)\\ab 

Here, it is most convenient to use directional orbitals along the a and b axes, with cL = — ^ c]^^ ± ^ cj^ . The 

density operators: hi(^ = cLc^^ and hi — hn^ + n^^, where (C|f } = 0, are restricted to single occupancies, 

and T^ ~ —\T^i i ^'^t ^^^° depend on the direction of the bond (ij). The pseudospin T == i orbital 
operators are given by 

Tl = \al = i(n„ - n,,), Tf = i< = i(4c,, + 4c„) . (5) 

Finally, the superexchange constant J = t^ jei^Ax) is given by the high-spin excitation energy e(M^i), see 
Ref. 1 3 1. The first term of Eq.Q favors FM spin and alternating orbital (AO) order, while the second one 
is optimized by AF spin and ferro orbital (FO) order — both have similar strength and therefore closely 
compete. 

The last term of the Hamiltonian is a crystal field splitting of Cg orbitals, 

H, = E,J2Tt = \e^ ^(n„ - n,,), (6) 

i i 

with the orbital splitting Ez- We consider only _E^ > with out-of -plane \z) orbitals favored, which is 
suggested by the elongated octahedra in LaSrMn04 along c axis [ 1 1. 

In order to treat this model, we employ a Markov chain Monte Carlo (MCMC) algorithm for the classical 
core spins (see, e.g. Ref. [6J) combined with exact diagonalization for the electron degrees of freedom 
|4][71. In this method, the many-particle Hamiltonian for each spin configuration S determined by {uij} 
is solved by exact diagonalization. The lowest states are then used to evaluate the trace over the fermionic 
degrees of freedom, Trc e^'^^''^' —: w{S), which gives the statistical weight for a given spin configuration 
iS, and which is sampled by the MCMC. Autocorrelation analysis was used in order to verify that enough 
configurations were skipped between measurements. The Boltzmann factors of the Lanczos eigenstates 
are measured in order to ensure that only negligible weight is lost when the weight is calculated from the 
lowest eigenstates and we found that at most one or two percent are missed for fit > 25, where /? = l/fc^T 
is inverse temperature. 

© 2003 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinheim 
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Fig. 1 Spin correlations along (0, l)-direction depending on distance r at various temperatures, as obtained with 
VS X \/8 clusters in ah plane for: (a) J' = 0, (b) J' = 0.02f, and (c) J' = 0.05f. Parameters: J/t = 1/8, and 




Fig. 2 Orbital correlation along (0, l)-direction depending on distance r at various temperatures, as obtained with 
\/8 X \/8 clusters in ah plane. Parameters and temperatures as Fig.Q 



3 Numerical results and discussion Spin correlations S(r) — {Si ■ Si+r) for the undoped system 
along the (0, l)-direction are gradually weakened when temperature increases (Fig.^. By increasing J' 
one identifies two distinctly different cases: for J' = the system is FM at low temperature f3t = 100 
[Fig. |l(a)|, wh ile AF spin correlations arise for J' — 0.02t and J' = 0.05i, leading to zig-zag lines in Figs. 
|l(b)| and |l(c)| Already a small value J' = 0.02t suffices to switch the spin order from FM to AF, because 
FM and AF terms strongly compete in Eq. @. For the same reason, increasing temperature can easily 
destroy the magnetic order and both FM and AF correlations decrease rapidly around pt ^ 50, except for 
the stronger AF superexchange J' — O.Q5t, with more robust AF coiTelations. 

Next we investigate orbital correlations T{r) = (TiTi^r) shown in Fig.|2l which exhibit complemen- 
tary behavior to the spin sector. Here we selected Ti = -j^{Tl + Tf) operators which correspond to the 

expected AO order with -7^{\x) ± \z)) orbital states in the 2D ab plane. One finds that the orbital correla- 
tions alternate for the FM spin order, while this alternation is suppressed for increasing J', as almost only 
orbitals within the plane are then occupied (FO order) which promotes the AF spin order. For J' = 0, see 
Fig. |2(a)| one notes that the AO order remains robust with rising temperature although the spin correlations 
vanish, see Fig. |l(a)| Remarkably, in the AF case one observes the gradual development of AO order as 
the AF spin correlations become weaker. For J' — Q.02t this crossover happens at fairly low temperature 
Pt — 50, while it moves to higher temperature for J' = 0.05t. 

As the magnetic order is so closely related to the orbital order in ab planes, one expects that changing 
orbital occupation will modify the magnetic order. Therefore, we investigated the electron density within 
the out-of-plane \z) orbitals, n^ — {niz). Indeed, n^ — 0.5 in the FM phase at low temperature (J' = 0), 



© 2003 WILEY-VCH Vcrlag GmbH & Co. KGaA, Weinheim 



1.0 



0.4 



0.2 



Daghofer, Neuber, Oles, and von der Linden: Temperature effects in a spin-orbital model for manganites 



(a) 



e- 



V- 



-^^J'=0 

- 0- J'=0.02 t 

-V J'=0.05t 



-o 



,J'=0.05t, E =0.2 t 



40 



60 



80 



100 



0.05 

0.04 

_< 0.03 

0.02 

0.01 
-0.01 

-0.02 

S-0.03 

-0.04 

-0.05 



(b) 




,::.--^ 


,^.-v---- 


G--" 


w- — r" 






-D 







(c) 










^.. 


:=V: 


:l::-- 




::_:,:.:g 


G- 





--Q--" 




- 













40 



60 



80 



100 



P 



Fig. 3 Temperature dependence of: (a) density in the out-of-plane \z) orbitals, as well as average eg superexchange 
terms: (b) Jap and (c) Jpm- Remaining parameters as Fig.|l(c)| 



consistent with AO order, while FO order with \x) orbitals (uz ^ 0) takes over for the strongly AF case 
(J' = 0.05t). With rising temperature, \z) orbitals are slightly depleted in the ferromagnet as orbital 
alternation becomes weaker. On the contrary, their occupation grows in the AF case, when such alternation 
develops. This happens quite soon for J' = 0.02t, but takes much longer for J' — 0.05t, where AF 
spin order also persists to higher temperature, see Fig. |2(c)| Further, if one adds a crystal field splitting 
Ez — 0.2t in the AF phase with J' — 0.05t, this forces almost all electrons into the out-of-plane \z) 
orbitals, the situation encountered in LaSrMn04. With rising temperature, however, electrons redistribute 
and in-plane \x) orbitals gradually fill in, see Fig. |3(a)| which was also observed in experiment Q]. 

The average AF and FM terms in eg superexchange can be derived from Ref. f^ and are given by: 
Jaf = ^J{nicnjc) - \J{'^T^TJ - ^n,hj), and JpM = ^J{2T^TJ - \hihj). Although Jfm becomes 
slightly smaller at higher temperature in the FM phase at J' = [Fig. HJb)], the changes are not very 
marked, because AO order persists. The FM coupling is therefore still almost the same as at low temper- 
ature, even if the magnetic order is lost because of thermal fluctuations. In the AF systems (J' = 0.02i, 
0.05t), the FM terms are first quenched by FO order, but become more important with rising temperature 
while the AF ones are reduced and, as before, this happens much sooner for J' — 0.02i than for J' = 0.05i. 
Interestingly, the AF eg couplings are then changed into FM ones with rising temperature [Fig.|3jc)]. 

Summarizing, at finite crystal field Ez ^ 0.2i, as in LaSrMn04, one observes that both JpM and Jaf 
are weaker than for Ez — 0. The reason is that pairs of alternating orbitals are missing for FO order with 
\z) orbitals in ah plane (small Jfm), and the overlap between \z) orbitals is much smaller than between \x) 
orbitals (small Jaf)- Therefore, the AF order observed in LaSrMn04 is promoted mainly by J'. 
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